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Abstract 

Background: With the migration of human beings, advances of agricultural sciences, evolution of planting patterns and 
global warming, soybeans have expanded to both tropical and high-latitude cold regions (HCRs). Unlike other regions, HCRs 
have much more significant and diverse photoperiods and temperature conditions over seasons or across latitudes, and 
HCR soybeans released there show rich diversity in maturity traits. However, HCR soybeans have not been as well classified 
into maturity groups (IVlGs) as other places. Therefore, it is necessary to identify MGs in HCRs and to genotype the maturity 
loci. 

Methods: Local varieties were collected from the northern part of Northeast China and the far-eastern region of Russia. 
Maturity group reference (MGR) soybeans of IVlGs MGOOO, MGOO, and MGO were used as references during field 
experiments. Both local varieties and IVIGR soybeans were planted for two years (2010-201 1) in Heihe (N SO'IS', E 127 '27', H 
168.5 m), China. The days to VE (emergence), Rl (beginning bloom) and R7 (beginning maturity) were recorded and 
statistically analyzed. Furthermore, some varieties were further genotyped at four molecularly-identified maturity loci El, E2, 
E3 and £4. 

Results: The HCR varieties were classified into MGO or even more early-maturing. In Heihe, some varieties matured much 
earlier than MGOOO, which is the most early-maturing known MG, and clustered into a separate group. We designated the 
group as MGOOOO, following the convention of MGs. HCR soybeans had relatively stable days to beginning bloom from 
emergence. The HCR varieties diversified into genotypes of El, E2, E3 and E4. These loci had different effects on maturity. 

Conclusion: HCRs diversify early-maturing MGs of soybean. MGOOOO, a new MG that matures much earlier than known MGs, 
was developed. HCR soybean breeding should focus more on shortening post-flowering reproductive growth. El, E2, E3, 
and E4 function differentially. 
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Introduction 

Soybean, a short-day crop of significant economic and 
ecological importance, diversifies significaiitiy during maturity. 
Maturity traits not only include total growth duration but also the 
vegetative and reproductive phases, and their relative ratios [1,2], 
reflecting whether temperature and light conditions satisfy the 
growth of soybean [3] and determine the adaptation of soybean 
varieties. 

Maturity is controlled by multiple loci or genes [4]. Nine 
maturity loci {E1-E8, and J) have been reported [5-13], have 
different effects on soybean flowering and play different roles in 
photo-thermal responses [9,14,15]. Four of these loci have been 



molecularly identified with map-based or candidate-based cloning. 
El encodes a soybean-specific potential transcription factor [16], 
E2 is a GIGANTEA homologxie of GmGIa [1 7], and E3 and E4 are 
the GmPhyA3 [18] and GmPhyA2 [19] phytochromes, respectively. 
More of these have been reviewed by Xia et al. [20] . Although 
some flowering genes are not included in or molecularly identified 
as maturity loci, they might also function in maturity. For 
example, GmFT2a displays differential transcriptional profiles 
under different temperatures and photoperiod conditions in two 
varieties with totally different photoperiod sensitivities [21]; 
however, its polymorphism appears to not be related to the 
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Table 1. Soybean varieties and their maturity groups (MGs) 



MGR (maturity group reference) Varieties from high-latitude cold regions with 

varieties proposed MG 



Maturity group North America China High-latitude cold region in Northeast China Amurskaya Oblast in Russia 



MGOOOO* Dongnong 36, Dongnong 41, Dongnong 41 -C*, Paula^, R-3'', R-4, Star 4/75 '', Sunset'' 

Hujiao 07-2123, Hujiao 07-2479'^, Lingbei s'' 



MGOOO 


Maple Presto, OAC 
Vision 


Heihe 12 


Dengke 2, Ha 6223-4, Heihe 11, Heihe 14, Heihe 20, 
Heihe 28, Heihe 35, Heihe 41, Heihe 44, Heihe 49, 
Kennong 8 


Bista, R-2, Sonata 


MGOO 


Canatto, Maple Ridge, 


Heihe 3, 


Bel 02-7495, Beidou 16, Beidou 24, Beifeng 1, 


Amur 262, Amur 283, Dewdrop, 




Glacier 


Heihe 8 


Dongnong 40, Fengshou 23, Hefeng 37, Heihe 5, 


Gritiaz 80, Harmony, Lydia, October 








Heihe 7, Heihe 9, Heihe 13, Heihe 33, Heihe 37, 


Revolution 70, R-1, Sunset 1, Terek 








Heihe 39, Heihe 45, Heihe 50, Heihe 51, Huajiang 2, 










Jiufeng 7, Jiufeng 10, Mengdou 7, Mengdou 9, Mengdou 11, 










Mengdou 31 




MGO 


MN0201, MN0901, 
Surge, Traill 


Jilin 30, 
Jiunong 21, 
Suinong 14 


Bel 1249, Beidou 8, Beidou 19, Beidou 37, Fengshou 15, 

Fengshou 24, Fengshou 26, Fengshou 27, Heihe 18, 

Heihe 36, Heihe 38, Heihe 43, Heihe 46, Heihe 48, Huajiang 3 





Jiangmodou 1, Jiufeng 9, Mengdou 30, Zhongzuo GHJ90962 



*MGOOOO is newly proposed here. 

^Proposed MG reference soybeans for MGOOOO. 

doi:l 0.1 371/journal.pone.00941 39.t001 



maturity diversification [22]. Therefore, the underlying mecha- 
nism of maturity diversity is not yet clear. 

Categorizing soybeans into difiFerent "maturity groups (MGs)" is 
convenient for breeding practices [23]. In North America, a 13- 
maturity group classification system has been set up according to 
the latitudes of adaptation [24—26]. In contrast, Chinese soybean 
researchers have divided varieties into twelve MGs based on the 
environments and planting patterns in China [27-29]. However, 
this classification method is not unified. Gai et al. classified 
Chinese soybeans into twelve maturity groups based on a maturity 
group distance of 10-15 days [29]. Wu et al. used the average of 



the neighboring maturity groups as the threshold to classily 
maturity groups [30], whereas Alliprandini et al. used a linear 
regression method to categorize local soybeans [31]. These 
different methods support the idea that the maturity mechanism 
is rather complex in soybeans and that diversified environments 
have significant impacts. However, in high-latitude cold regions 
(HCRs), littie is known. 

Since the development of soybean breeding and plantations, 
HCRs have become important for soybean production. Soybean 
production HCRs are predominantly located in China, Russia and 
North America. In China, they are located in the HeUie, Yichun, 



Table 2. Days to beginning maturity from emergence for the MGR (maturity group reference) soybeans in Heihe, China. 









Days to begin 


ning maturity from emergence 


Variety 


MG 


2010 


2011 


Maple Presto (E) 


MGOOO 


88.5 ±0.9 


82.0 ±1.4 


OAC Vision (L) 




93.8±0.4 


87.8±0.9 


Heihe 12 




94.6 ±0.9 


85.9±1.5 


Canatto (E) 


MGOO 


97.8 ±0.9 


92.0±0.8 


Maple Ridge (E) 




95.9±0.8 


96.3 ±0.6 


Glacier (L) 




107.3±1.1 


95.2±0.8 


Heihe 3 




102.4±1.2 


97.8±0.9 


Heihe 8 




97.3±1.1 


90.8±0.6 


MN0201 (E) 


MGO 


120.4±1.0 


99.4±0.7 


Traill (E) 




1 1 3.1 ± 1 .2 


99.1 ±0.7 


MN0901 (L) 




124.9 ±1.2 


115.4±0.6 


Surge (L) 




R6 


R5 


Suinong 14 




R6 


R5 


Jllln 30 




R4 




Jiunong 21 



Soybeans with brackets are maturity group reference soybeans from North America. The others are from China. Letters E and L (in brackets) indicate that the soybean is 
relatively early- or late-maturing, respectively, in the corresponding maturity group; R4, R5 and R6, the growth stages of soybeans failing to mature at first frost. 
doi:l 0.1 371/journal.pone.00941 39.t002 
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Figure 1. Days to beginning maturity and beginning bloom from emergence in 2010. (a). Days to beginning maturity from emergence 
(2010). (b). Days to beginning bloom from emergence (2010). MGOOOO, a new proposed maturity group, is marked in red. 
doi:1 0.1 371 /journal.pone.00941 39.g001 



Daxing'anlmg region, North Qjqiha'er and North Suihua of 
Heilongjiang Province and Hulunbei'er of Inner Mongolia, which 
are predominantly located between N47°— 53°34'. In these 
regions, soybean acreage reaches over 1,000,000 ha [32]. In 
Russia, the far-eastern region, especially the Amurskaya Oblast, 
which neighbors the Heilongjiang province of China, is the main 
region for soybean production [33], accounting for 90% of 
Russian soybean production [34]. In Canada, soybean production 
has expanded from southern Ontario in the 1970s to eastern 
Ontario, Quebec and Manitoba, reaching 1.7 million ha in 
2012 (http://www.soybeancouncLl.ca/). HCRs will become more 



important for soybean production, especially under the back- 
ground of climate change. 

In HCRs, the environment changes more violently across 
latitudes and between seasons than at lower latitudes, especially for 
temperature and photoperiod. In these areas, soybean has a short 
growth duration and is relatively insensitive to photo-thermal 
regimes [35]. Here, we use MGR soybeans, MGO, MGOO and 
MGOOO, as a reference to classily local HCR soybeans into 
diflferent maturity groups and to genotype them at four 
molecularly identified maturity loci. El, E2, E3 and E4, to learn 
the adaptation mechanism of soybean in extreme environments. 



Table 3. Soybean varieties earlier than the MGR (maturity group reference) soybeans of MGOOO in terms of days to beginning 
maturity from emergence. 





Variety 


2010 


2011 


Variation/Mean 


R-4 


78.3±0.5 


72.1 ±0.9 


0.08 


Star 4/75* 


78.2±0.5 


73.3±1.3 


0.06 


Hujiao 07-2123 


79.1 ±1.2 


74.0± 1 .4 


0.07 


Sunset* 


77.9±0.8 


742±0.9 


0.05 


Hujiao 07-2479* 


78.9±1.0 


76.3±0.5 


0.03 


R-3* 


80.0±0.5 


76.6±0.9 


0.04 


Dongnong 36 


85.5±0.9 


77.3±0.7 


0.10 


Paula* 


79.8±0.7 


77.8±1.6 


0.03 


Dongnong 41 


84.3±1.1 


77.9±0.6 


0.08 


LIngbel 8* 


82.5±0.7 


79.1 ±0.7 


0.04 


Dongnong 41 -C* 


84.0±0.8 


81.8±0.7 


0.03 


MGR soybeans, MGOOO 


88.5-94.5 


81.9-87.8 


NA 



*MGR soybeans for proposed MGOOOO. NA, "not available". 
doi:l 0.1 371/journal.pone.00941 39.t003 
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Materials and Methods 

1 Plant materials 

Ninety-three soybean varieties were involved in tliis study 
(Table 1). Of these varieties, 9 were MGR soybeans (MGOOO, 
MGOO, and MGO) from North America; 6 were MGR varieties 
(MGOOO, MGOO and MGO) from China; 18 were from Amurskaya 
Oblast (Russia); and the remainmg 60 were released varieties of 
HCRs from China. 

2 Field experiments 

Field experiments were conducted between 2010 and 2011 in 
the experimental field of Heihe Branch of HeUongjiang Academy 
of Agricultural Sciences (N 50°15', E 127°27', H 168.5 m). No 
specific permissions were required for these locations/activities. 
No endangered or protected species were involved. On May 9* of 
each year of these two years, the soybeans were manually sowed in 
two rows (row length of 2 m, plant spacing of 5 cm, and row 
spacing of 60 cm). After emergence, they were thinned for 30 
uniform healthy plants. The days reaching VE (emergence), Rl 
(beginning bloom) and R7 (beginning maturity) were recorded in 
line with Fehr and Caviness [36]. 

3 Genotyping of El, E2, E3 and E4 

Fresh leaves were used to isolate DNA using the TIANGEN 
(Beijing, China) New Plant Genomic DNA Isolation kit (DP320- 
03). Loci El, E2, E3 and E4 were genotyped following previously 
described methods [16-19]. 

4 Statistical analysis 

In R environment [37], APCluster [38] was used to cluster 
varieties into different maturity groups with negative Euclidean 
distances as mutual pairwise similarities and with the maturity 
duration data in 2010 and 201 1 as input. A principal component 
analysis (PCA) was also conducted in FactoMineR [39] with 
default settings. 

Results 

Tested soybean varieties matured differently with 
different beginning bloom days 

A high latitude with a long photoperiod and low average 
temperature significantly affected the maturations rates of the 
soybeans (shown in Table SI and File SI). For the MGR soybeans 
(MGOOO - MGO), MGO diverged with some varieties maturing 
and some not, while MGOO and MGOOO matured (Table 2). The 
matured MGR soybeans matured 82.0 to 124.9 days after 
emergence (Table 2), and their MGs significantly differed at the 
day to maturity after emergence (Figure la). Moreover, the days to 
beginning bloom from emergence (DFF) also diversified signifi- 
candy (Figure lb). MG was the more early-maturing with the 
range of DFF narrowed even more (Figure lb). 

For the other non-MGR soybeans, they also exhibited 
significant variations in their maturity, with a range of 43-54 
days (Table SI). Huajiang 3 (2010), and Beidou 37 and Zhongzuo 
GHJ90962 (2011) matured last in different years, whereas Sunset 
(2010), and R-4 (2011) were the earliest to mature (Table SI). 
However, 1 1 varieties matured even ealier than MGOOO during 
the two years, that is, R-4, Star 4/75, Hujiao 07-2123, Sunset, 
Hujiao 07-2479, R-3, Dongnong 36, Paula, Dongnong 41, Lingbei 
8, and Dongnong 41-C (Table 3). Dengke 2 matured earlier than 
the MGR soybeans MGOOO in 2010 but later than the early 
MGOOO reference varieties in 2011 (Table SI). 
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Figure 2. Heatmap for the clustering of soybeans from high-latitude cold regions. 

doi:10.1371/journal.pone.0094139.g002 



The tested soybean varieties clustered into four maturity 
groups 

The tested soybean varieties consistently matured between years 
(Table 4). R1-R7 and VE-R7 displayed a high correlation in the 
same year, whereas VE-Rl did not correlate with the former two 
(Table 4). With APCluster [38], these soybeans were clustered into 
seven clades (Figures 2 and 3). For the MGR soybeans, MGO and 
MGOO each clustered into two clades, whereas MGOOO clustered 
into one clade (Figures 2 and 3). The varieties that matured earlier 
than MGOOO during the two years all formed one clade (Figures 2 
and 3). This clade is significantly different with the other MGs and 
has been designated as a new MG of MGOOOO following the 



convention of MG. The other clade, without MGR soybeans, is 
located between the two MGO clades; thus, it is also included in 
MGO. Thus, the tested varieties were clustered into four MGs 
(MGO, MGOO, MGOOO, and MGOOOO, see Table 1). Based on the 
Variation/Mean (VE-R7) and the VE-R7 during the two years, 
MGR soybeans for MGOOOO were proposed. They are Star 4/75, 
Sunset, Hujiao 07-2479, R-3, Paula, Lingbei 8, and Dongnong 41- 
C (Table 1). 

Genotyping of the El, E2, E3 and E4 maturity loci 

Forty-eight varieties were genotyped for the maturity loci El, 
E2, E3 and E4 (Table 5). Of those varieties, 7, 8, 21, and 12 were 
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Figure 3. Clusters of soybeans from high-latitude cold regions. Numbers from 1 to 10 and 12 indicate different maturity group reference 
soybeans. 
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from MGOOOO, MGOOO, MGOO, and MOO, respectively. Seven 
genotypes were identified, that is, ele2e3e4 (4 varieties in 3 MGs), 
eIe2e3E4 (23 in 3), ele2E3e4 (4 in 2), ele2E3E4 (12 in 3), 
elE2E3E4 (1 in 1), Ele2e3e4 (2 in 1), and Ele2e3E4 (2 in 1). Each 
MG diversified into El, E2, E3 and E4 genotypes. 

El, E2, E3 and E4 maturity loci are diverse in maturity 
groups 

PCA showed that the first two principal components of Dim 1 
and 2 could explain 71.27% and 13.60% of the variation, 
respectively, and 84.87% in total (Figures 4 and 5). During the two 
years, VE-R7, VE-Rl, and R1-R7 all showed positive correlations 
with the frrst principal component. Dim 1, whereas VE-R7 had 
the highest correlation (Figure 4). As for the second component. 
Dim 2, VE-Rl positively correlated with it while R1-R7 
negatively correlated (Figure 4). That is to say, Dim 1 reflects 
maturity while Dim 2 relates with flowering. In Figure 5, the 
maturity groups MGOOOO, MGOOO, MGOO and MGO dispersed 
along Dim 1 . On Dim 1 , the maturity loci El and E4 had larger 
coordinates than their recessive versions, el and e4, whereas E2 
and E4 were smaller than their recessive versions, e2 and e4. 

Discussion 

Until recently, HCRs were not viewed as suitable for soybean 
production because they exhibit low average temperatures and 
short frost-free periods that are not long enough for late-maturing 
varieties. However, with the variety improvement and develop- 
ment of farming technology, HCRs have become more and more 
important for the soybean industry under the background of global 
warming, the increasing population, and arable land shortages. 

The soybean is a short-day crop and diversifies, in terms of 
maturity, among varieties. For convenience, soybean varieties are 
classified into different MGs to direct breeding and planting 
practices. However, although HCR soybeans are known to be 
relatively insensitive to photoperiods compared with those from 



other regions, information is stiU missing regarding their MGs. 
This information is missing partially because HCR soybeans 
mature much earlier and have a low differentiation degree at 
medium and low latitudes, where the main production areas 
currently are. However, as we know, few reports about MG have 
been conducted in HCRs. 

In Heihe, the tested varieties displayed a broad diversity in their 
maturity (Figure 1, and Tables 2 and SI). MGO matured partially, 
suggesting that HCR soybeans belong to MGO or even more early- 
maturing MGs. Some HCR soybeans matured much earlier than 
MGOOO and were clustered into the new MGOOOO. Therefore, 
HCR soybeans show diverse maturity rates. 

HCR soybeans are relatively insensitive to the photoperiod. 
Thus, whether they can flower is not the main problem for 
maturation. In Table 4, the correlation between VE-Rl and VE- 
R7 or R1-R7 was less than that between VE-R7 and R1-R7, 
suggesting that the reproductive growth duration is much more 
important than the vegetative growth/duration for maturation. 
The PCA results also provide further evidence. The first principal 
component, Dim 1, related most significantly with maturity, 
whereas the second component. Dim 2, exhibited flowering time 
(Figure 4). The tested varieties predominandy dispersed along Dim 
1 , whereas they narrowed along Dim 2 (Figure 5), suggesting that 
the post-flowering period is much more important than the pre- 
flowering period in maturity diversification and classification and 
that post-flowering photoperiodic responses play an important role 
for HCR soybean maturation [40]. 

Soybean growth period traits not only include the maturity time 
but also its structure [1], i.e., the time to first-flowering, first- 
podding, etc., should also be considered. Egli found that seed yield 
is related to the length of the reproductive phase rather than the 
total growth duration [2]. For early maturing varieties, SD 
treatment before flowering significantly promotes post-flowering 
development but does not significantiy promote pre-flowering 
development [41]. Our experiment also yielded similar results. In 
our experiment, Paula and Ha 6223-4 had the top two shortest 
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Table 5. Genotype of soybean varieties at the El, E2, E3 and E4 loci. 





Maturity Group 


Genotype 


Number 


Variety 


MGOOOO 


ele2e3e4 


2 


Hujiao 07-2123, Hujiao 07-2479 




ele2E3e4 


2 


Dongnong 36, Paula* 




ele2E3E4 


2 


R-4*, Star 4/75* 




elE2E3E4 


1 


Sunset* 


MGOOO 


e1e2e3e4 


1 


Dengke 2 




ele2e3E4 


1 


Heihe 12 




ele2E3e4 


2 


Bista*, R-2* 




ele2E3E4 


4 


Heihe 14, Heihe 35, Heihe 49, Sonata* 


MGOO 


e1e2e3e4 


1 


Mengdou 31 




ele2e3E4 


12 


Amur 262*, Dongnong 40, Fengshou 23, Hefeng 37, Heihe 3, Heihe 8, Heihe 9, Heihe 33, Heihe 45, 
Huajiang 2, Mengdou 11, R-1* 




ele2E3E4 


6 


Amur 283*, Dewdrop*, Heihe 13, Lydia*, October Revolution 70*, Terel<* 




E1e2e3e4 


2 


Harmony*, Gritiaz 80* 


MGO 


ele2e3E4 


10 


Beidou 37, Fengshou 15, Fengshou 26, Fengshou 27, Heihe 18, Heihe 38, Heihe 43, Huajiang 3, 
Mengdou 30, Sunset 1* 




E1e2e3E4 


2 


Jiangmodou 1, Heihe 36 



*Russian varieties. 

doi:l 0.1 371 /journal.pone.00941 39.t005 



that HCR soybean breeding should guarantee sufficient vegetative 
growth. Under high latitudes with low average temperatures and 
short frost-free periods, soybean breeding should be more focused 
on shortening post-flowering reproductive growth when vegetative 
growth is shortened to an extent that necessary vegetative 
accumulation can be guaranteed. At low latitudes or in tropical 
and subtropical settings, a short photoperiod is too promotional for 
soybean flowering and maturation, and a long juvenile trait is 
needed to lengthen the vegetative growth period to ensure enough 
vegetative accumulation to increase production [42,43]. There- 
fore, in breeding, it is necessary to regulate vegetative and 
reproductive growth to adapt to different environments. 

Some varieties matured much earlier than the MGR soybeans 
of MGOOO, which is known to be the earliest maturing MG 
(Table 3 and Figures 2 and 3). The distance between these early 
varieties and the MGR soybeans of MGOOO is sufficient to 
determine a new MG, according to the conventions of MG 
classification. Star 4/75, Sunset, Dongnong 36, and Dongnong 41- 
C should be the MGR soybeans of this new MG because they 
have a relatively low SD/Mean value (Table SI). 

Soybean maturity traits are under the control of multiple genes. 
Currently, nine related loci have been found [5-13], and their El, 
E2, E3 and E4 have been molecularly identified [16-19]. Their 
reproductive periods were demonstrated to be controlled by one 
major gene, plus polygenes, in a research study involving nine 
crosses, and two major effector QTLs, qRP-c-1 and qRP-l-1, (which 
associate with E8 and E3, respectively) may function there [44] . In 
the 48 tested varieties, loci El and E2 were found to exist 
predominantly as their recessive alleles, and their dominant alleles. 
El and E2, only distributed in 4 and 1 varieties, respectively 
(Table 5), suggesting that loci El and E2 have strong impacts on 
the adaption to HCRs. As for loci E3 and E4, they appear to have 
different roles on maturity. In Figure 5, allele E3 has a lower 
coordinate on Dim 1 than its recessive allele, e3, whereas allele E4 
has a higher coordinate than its recessive allele, e4. Thus, allele E3 
promotes maturation while allele E4 delays it. From their 
coordinates on Dim 2, it appears that E4/e4 does not have a 
role in flowering, whereas E3 delays and e3 promotes flowering. 



VE-Rl in 2010 (Table SI), but belonged to different MGs 
(MGOOOO and MGOOO in Table 1, respectively), and MGOOOO had 
the narrowest range of DFF (Figure lb). These observations 
suggest that HCR soybeans have relatively stable DFF. Therefore, 
for HCR soybeans, which are typically early maturing, their 
vegetative growth periods are similar, but their reproductive 
periods are diverse (Table SI). This further suggests that 
lengthening the reproductive growth period would not hamper 
the vegetative growth nor inhibit dry matter accumulation [1] and 



Variables factor map (PCA) 




-0.5 0.0 0.5 

Dim 1 (71.27%) 



Figure 4. PCA variable factor map of soybeans from the high- 
latitude cold regions. 2010. VE-R7 is the days between VE and R7 in 
2010. Others are similar. 
doi:1 0.1 371/journal.pone.00941 39.g004 
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Individuals factor map (PCA) 
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Figure 5. PCA individuals' factor map of soybeans from high-latitude cold regions. E1/e1, E2/e2, E3/e3, and E4/e4 are the four maturity 
loci marked in red. MGO, MGOO, MGOOO and MGOOOO are the four maturity groups marked in red. The top ten contributor soybeans are marked in 
black. 

doi:1 0.1 371 /journal.pone.00941 39.g005 



The converse roles of two phytochrome genes, E3 and E4, might 
be related with the low average temperatures in HCRs. In 
addition, each MG had multiple genotypes, and one genotype 
could distribute into several MGs. Thus, the number of recessive 
alleles did not correlate with the MG classification in the HCRs. 
These four MGs harbor the diversity of genotypes of E genes, 
indicating that, in HCRs, the mechanism of maturity is rather 
complex and involves more genes than these four E genes. 

Supporting Information 

File SI Raw data of soybean maturity traits in 2010 and 201 1. 
(XLS) 

Table SI Growth periods of the tested soybean varieties. 
(DOC) 
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